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q ■ Abstract 



We prove the existence of ground state in a multidimensional nonlinear Schrodinger model 
of paraxial beam propagation in isotropic local media with saturable nonlinearity. Such 
ground states exist in the form of bright counterpropagating solitons. From the proof, 



03 



a general threshold condition on the beam coupling constant for the existence of such 
fundamental solitons follows. 



1 Introduction 



> 

^ ■ The existence of solitary waves in nonlinear evolution partial differential equations has been 
CNj ■ the major concern since the beginnings of the field 0. The existence of ground state 
^ . or fundamental soliton in one-dimensional nonlinear Schrodinger (NLS) equation was never 
0^ '. much of a concern - their existence and stability followed from the inverse scattering theory 
jl]. However, in more than one dimension this was, and still is, an open question ||. In 
^ ■ multidimensional cases there exists no mathematically rigorous theory that would guarantee 
J> ■ their existence and uniqueness - let alone stability in propagation. 

In this paper we prove the existence of ground state in the form of counterpropagating 
^! solitons in an isotropic local saturable NLS model of beam propagation |4j]. This model rep- 
resents a physically relevant description of soliton generation, following from the theory of 
photorefractive effect in crystals that respond to light by changing their index of refraction. 
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2 The model 



We consider the following dimensionless NLS system: 



iF z + A ± F + TE F = 0, 
-iB z + A ± B + TE B = 
d t E + E { 



io (2.1; 

J = I^VlF! 2 , 



where F and i? are the slowly-varying envelopes of the forward and backward propagating 
beams, z (0 < z < L) is the propagation coordinate, a; = (xi,X2) € IR 2 are the transverse 

2 

coordinates, and Aj_ = is the transverse Laplacian. Furthermore, V is the beam coupling 

constant, F the homogenous part of the space charge field generated in the photorefractive 
crystal, and Iq is the beam intensity, expressed in terms of the background intensity The 



j I F III B I 

time independent solution of (|2.1|) must satisfy d t E = 0, i.e. F = — = — ■ We 



consider such a situation. Then (12.11) becomes 



^ + A^F - T^j^F = 0, 
-iB z + A ± B - r ^+^ B = o. 

Here, we assume the following boundary and initial conditions. 
Boundary conditions: F (x, z) , B (x, z) — > as \x\ — >• oo, for z < L. 
Initial conditions: F (x, 0) = F (x) , 5 (x, L) = 5l (x) are given. 
Equations ( |2 .2| ) can be written as iF z = 5E ^ F ^ , — zF z = 5E \^ F ^ , where 

F[£,F] = / |(|V J B| 2 +|VF| 2 ) + ir[| J B| 2 + |F| 2 -ln(l + |F| 2 + |F| 2 )] 

is the energy functional. The integral (here and elsewhere) is taken across the whole transverse 
plane. One basic conservation law of the system (|2.2|) is the power conservation, given by 



(2.2) 



|£| 2 + |F| 2 = / |F | 2 + |F | 2 for < z < L. (2.3) 



To get the standing wave profiles of the system fl2.2f) , we set 

F (x, z) = e iXz u (x) , B (x, z) = e~ iXz v (x) , 
and then the system fl2.2|) can be transformed into the following system: 



,2i„2 



A±u - r " tl 2 u = Xu, 
- T ,?+? 2 v = Xv, 

so that the condition (0.31) also can be normalized as follows: 



(2.4) 



P[u, v] — / (y 2 + f 2 )=l (normalization). 



(2.5) 



Here, A is the propagation constant, which can be regarded as the chemical potential in other 
physical settings; mathematically it is the Lagrange multiplier following from condition ( [2.51 ). 
The system of equations ( |2.4j ) has been solved numerically in [EJ, but with an external lattice 
potential I g included instead of the uniform background intensity = 1. Fundamental coun- 
terpropagating solitons have been obtained and a threshold condition determined. However, 
these results have been numerical, without rigorous proofs to substantiate their existence. 

To obtain the ground state of ( |2.2|) rigorously, we consider the following energy minimization 
problem: 

fir = inf {E [u, v] : u,v G H 1 (IR 2 ) , P[u,v] = l}, (2.6) 
where the energy functional now is 

E [u, v] = [ (| Vu\ 2 + \Vv\ 2 ) +T [u 2 + v 2 - In (l + u 2 + v 2 )] . 
Jr 2 



The main issue with the problem ( |2.6|) is whether the value pr (depending on the coupling 
constant T) can be achieved at a minimizer called the ground state solution of ( |2.2p . It is obvious 
that E [u, v] > — |T| for u, v G H 1 (M 2 ) , / (u, v) = 1, which implies the value fir > — |T| > — oo. 
Note that for the NLS equations with power nonlinearity, the infimum energy may not exist 
for some power magnitudes (see [0). Here, as T < 0, i.e., the self-focusing case, the potential 
energy of the saturable nonlinearity 



/ r> 2 -/n(i + P 2 )] 

Jr 2 



may compete with the kinetic energy 



|Vp| 
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so the magnitude of T affects the existence of the ground state solution. The main results may 
be stated as follows: 

Theorem 2.1. Let To be the following positive constant: 

T = inf =fe |VW|2 

weH 1 (IR 2 ) Jr 2 ^ 2 ~ ln ( 1 + w2 )} 

IMI 2 = 1 

(%) If V > —Tq, then fir can not be attained by a minimizer i.e. there is no ground state 
solution. 

(ii) If V < —Tq, then fir < and there exists a ground state solution which is radially 
symmetric and is denoted by (u,v) = p (r) (cos (p, sin (ft), where <p G IR is an arbitrary 
constant and p = p(r) is the energy minimizer of the following problem: 

Minimize H [p] over p G H 1 (R 2 ) , / p 2 = 1, (2.7) 

Jr 2 

where 

H[ P }= I \Vp\ 2 + T[p 2 -ln{l+p 2 )}. (2.8) 
Jr 2 
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The physical meaning is as follows: Theorem [2.1| indicates that ground states only can 
behave like bright solitary waves in saturable photorefractive media. The constant — T is the 
threshold for the existence of ground state solutions, which may be changed under the effect of 
external intensity lb (see Theorem [2.1| (ii)). 



Remark 2.2. By Schwartz symmetrization, it is obvious that the minimizer p of the problem 
( $■4 ) must be radially symmetric and its Euler- Lagrange equation may be expressed as follows: 

'p" + }p'-T T ^ = \pforr>0, 
p'(0) = 0, p(0)>0, 



where A is the Lagrange multiplier of the problem ( \2.^j ). We shall show that A > and the 
minimizer p is a positive and monotone decreasing function which decays to zero exponentially 
as the variable r goes to infinity. 

3 Proof of Theorem gj] 

The proof is divided in a number of steps. The following lemma is crucial in proving Theorem 
Lemma 3.1. The value pr defined in ( \2.6\ ) satisfies 



p T = inf{H[p\: peH'iR 2 ), [ 



P 2 = l 



Proof. Let u = p cos cj) and v = p sin </>, where both p and are H 1 functions. Then 

\Vu\ 2 + \Vv\ 2 = |Vp| 2 + p 2 |V0| 2 
and u 2 + v 2 = p 2 , and hence the proof of Lemma |3.1| is obvious. □ 



3.1 Proof of Theorem gT| (i) 

Proposition 3.2. Let 
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Then T n > 0. 



1 1 tt7 1 ' 



To prove Proposition 3.2, we need 



Claim Al. sup *-M 1 ~ = I. 

s>0 

Proof. Let h (s) = 5 ~^ ( 2 1+s) for s > 0. Then h! (s) = -^rfe + %ln (1 + s) and {s 3 h' (s))' = 

— ,^ s \2 < for s > 0, which implies h! (s) < for s > 0, i.e. h is a monotone decreasing 

function for s > 0. Note that s 3 h' (s) = at s = 0. On the other hand, by direct calculation, 
Urn his) = h and hence the proof is complete . □ 



By Claim Al, w 2 — ln(l+w 2 ) < ^w A , so J R2 [w 2 — In (1 + w 2 )} < \ J R2 w 4 for w £ 



2 

H 1 (R 2 ) , || w || 2 = 1, which implies 



T > inf 

weHUR 2 ) 2IMI4 



\w\\ 2 = 1 



On the other hand, by (2.2.5) in ||, 

L 2 |Vt«| 2 
inf T „ 4 = 25 2 , 4 > , 

Hl 2 = i 

where <S , 2 4 is the Sobolev constant. Therefore, T > 2^4 > and we complete the proof of 



Proposition 3.2 



Proposition 3.3. Suppose T £ (— T ,0) z.e. > T > — T . T/ien f/ie wa/ue /ir can not attain 
a minimizer such that /ir < 0. 

Proof. We may prove by contradiction. Suppose there exists w a minimizer of the value fir 
such that /ir < 0. Then J Rd \Vu\ 2 + T J Rd [u 2 — In (1 + w 2 )] < and ||it|| 2 = 1. Hence 



/ |W| 2 <-r/ [ M 2 -/n(l + W 2 )]=-^(T / 









\Vu\ 


J I J I 



which implies u = 0, since < -jr < 1. However, w = contradicts ||w|| 2 — 1- Therefore, 



we 



have completed the proof. □ 

Proposition 3.4. Suppose T £ (— Tq,0), i.e. > T > —To. T/ien i/ie value /ir can noi attain 
a minimizer such that fir > 0. 

Proof. It can be proved by contradiction. Suppose there exists a minimizer tt of the value /ir 
such that /ir > 0. Then u satisfies the Euler-Lagrange equation of the problem 



/ir = inf{tf[p]: p £ H 1 (R 2 ) , ^ p 2 = lj 



given by 



u 3 



Aw-T - = Awforx£R 2 (3.1) 

1 + u 2 

with ||u|| 2 = 1 and u(x) — > as |x| — >■ 00, where A is the Lagrange multiplier. Multiplying 
equation (|3.1|) by x- Vu and integrating over R 2 , we may derive the Pohozaev identity as follows: 



X = -T [ [u 2 -ln{l + u 2 )] . 

</R2 



(3.2) 
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The derivation is quite standard, so we omit the details here. On the other hand, multiplying 
equation ( |3.1| ) by u and integrating over IR 2 , gives 



A=- l \Vu\ 2 -T I _^L. (3.3) 

JR 2 Jr 2 l + u 2 

Here we have used integration by parts. Suppose /ir > 0. Then it is obvious that 



r f [u 2 - In (l + u 2 )] > - I \Vu 
Jr 2 Jr 2 

Hence (|3.3| ) implies 

A<T / [u 2 -ln(l + u 2 )] -r / - 
Jr 2 Jr 2 1 



2 



+ u 2 



Combining with ( [3.2|) and (|3.4j ), we have 

-2r f [u 2 -ln(l + u 2 )] <-Y ! - 
Jr 2 Jr 2 1 



u 4 



+ u 2 



which is equivalent to 



since T < 0. Let 



Then F(0) = and 



s 2 



F (s) = 2 \s - In (1 + s)l for s > . 

w L v n 1 + s 



s 2 



F' (s) = ? > for s > , 

KJ (1 + s) 2 



and hence F (s) > for s > 0, i.e. 



s 2 



(3.4) 



2 /_> 2 -/n(l + . 2 )]< / ^ (3.5) 



2\s-ln(l + s)] > for s > . 

Therefore, replacing s by -u 2 , 2 J R2 [ti 2 — Zn (1 + w 2 )] > J R2 which contradicts ( p.5|) and we 
have completed the proof. □ 

Proposition 3.5. Suppose T > 0. Then a? = can noi attain a minimizer. 

Proof. Let w E H 1 (R 2 ) and ||u;|| 2 = 1. For 5 > 0, let w<j (x) = (5x) for x E IR 2 . 
Then Hu^l^ = HHI2 = 1 an d Jr2 |Vu\5 (x)| 2 = o" 2 L 2 I Vw (w)| 2 dw. Moreover, by Claim Al, 
/ K 2 [wf — ln(l + w 2 )] < \ J R2 wj = \5 2 J R2 (y) dy. Hence H [ws] = O (5 2 ) tends to zero as 
8 goes to zero. This implies that n-p = 0. On the other hand, since T > 0, it is obvious that 
uy = can not attain a minimizer. Therefore, the proof is completed. □ 



Combining Propositions |3.2| - [3.5| , completes the proof of Theorem |2.1| (i 
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3.2 Proof of Theorem [27T] (ii) 

For the proof of Theorem [2.1| (ii), we firstly consider the following problem: 

fi T> e = inf H £ [u] , 
u e Hi (b± 

Pe M = 1 

where 



H e [u]= I \Vu\ 2 + T[u 2 -ln{l + u 2 )] 



and P £ [u] = f Bi u 2 for e > and -u G Hq f-^ 1 )- Hereafter, Si is the ball with radius ^ and 
center at the origin. 

Lemma 3.6. Assume T < —T < 0. Then 

(i) For e > 0, /ir, e can be achieved by a minimizer u e = u £ (r) > with radial symmetry. 

(ii) For e > sufficiently small, «r,£ < — cq, where cq is a positive constant independent of e. 

Proof. Fix e > arbitrarily. Since < s 2 — In (1 + s 2 ) < Cs 2 for set, where C is a positive 
constant independent of s, then H £ [u] > YC f g ^ u 2 = TCP £ [u] = YC for u G Hi (bi J and 

£ 

P £ [u] = 1. Hence the value /ir, e = in/ H £ [u] exists. Let {vk}^ =1 be a minimizing 

u G Hi (b^J 

Pe [U] = 1 

sequence of the value /xr, e - Without loss of generality, we may assume [ufc] I /ir, £ as A; — > oo. 
Note that each G f ^-J anc ^ ^ ^ = ^' ^PP^ symmetry rearrangement on each 

v k , say v* k G Sq 1 ^B^J and P £ = P £ [v k ] = 1. Note that / Bi |V^| 2 < / Bj |Vv fe | 2 and 

I Bx _ M) 2 - In (1 + K) 2 )] = / fli [(t; fc ) 2 - Zn (l + (v k ) 2 )]. Then ^ < # £ [t£] < # e N for 

£ £ 

all fc. Hence we may replace i>fc by v% and regard f^'s as functions with radial symmetry. On 
the other hand, since H £ [v k ] i Ur, e a s k — > oo and < s 2 — in (1 + s 2 ) < Cs 2 for s G M, then 

/ \Vv k \ 2 < C u 

£ 

where C\ is a positive constant independent of k. Hence Sobolev embedding gives v k — > u £ 
weakly in Hi f^lj as A; — oo (up to a subsequence). Moreover, by the Sobolev compact 

embedding, v k — > u £ in L 2 I Bi I as k — > oo. Consequently, P £ [u £ ] = lim P £ [v k ] = 1 and 

V £ / fc ^°° 



iUr,£ < H £ [u e ] < lim inf H £ [v^] = Hr, £ - This implies that u £ is the minimizer of the value 

k— )-oo 

/Ur i£ . Here, we have used Fatou's Lemma. Moreover, since each is radially symmetric and 
Vk — > u £ in L 2 I B\ I as fc — > oo, then u £ is radially symmetric. This completes the proof of 



Lemma [3.6| (i). 

The proof of Lemma |3.6| (ii) is given as follows: From the definition of To, 

V<5 > 0, 3w 5 e H 1 (R 2 ) , \\w s \L = 1 such that T > -= { R2 ^ - 5 . 

J R2 w 2 5 -ln{l + w 2 ) 

Let 5 = -E^. Then 5 > and T < — T n — 5, i.e. -r > T n + 5 due to T < -Tn. Hence 



-r > . J f |V ?' ay , i.e. < 0. Let w Se = where <^> £ e (R 2 ) is a cut-off 

function such that: ip £ = 1 in £?i , <y2 £ = in B\ = R 2 — 5i and ||y3 £ || c i = ^ where O (1) is 

e~ 1 e 

a bounded quantity. Then ws >e £ Hq (^l) an d ||w<5,e|| 2 = 1- Note that e and 5 are independent 
of each other. It is obvious that J BC w 2 + \ Vwsf = o £ (1), where B\ = R 2 — Bi and o £ (1) 

£ 

is a small quantity tending to zero as e goes to zero. Moreover, Hu^^]^ = ||w,s|| 2 + o £ (1) = 1 + 
o £ (1), f R 2 |V (ww e )| 2 = j R2 \Vw s \ 2 (f 2 +j R2 2 (w 5 V«; 5 )-(^ e V^ e )+w||V^ e | 2 = j R2 \Vw s \ 2 +o E (1), 
and J R2 {(w S (p £ ) 2 - In [l + («We) 2 ] } = / R a [^ 2 -ln(l + w 2 )] + o e (1). Hence from H [w s ] < 0, 
we get H e [ws <e ] = H [w$] + o £ (1) < — c as £ gets sufficiently small, where c = —\H [wg] > 
is a constant independent of e. Therefore, /ir, £ < ^ e [wg,e] < — c o an d we complete the proof of 
Lemma B.6I (ii). □ 



Remark 3.7. From Lemma \3.(\ (ii), we get fir < £tr, £ < ~cq < 0. 



Lemma 3.8. Under the same hypothesis of Lemma\TQ, the minimizer u £ satisfies \\u £ \\ / \ < 

K /or £ > sufficiently small, where K is a positive constant independent of e. 
Proof. By Lemma [3.6| (ii) /ir, £ < — Cq, which implies 

-c > / |Vu e | +r I u 2 - In (l + u 2 £ ) > [ \Vu £ \ +r. 

£ £ £ 

Here we have used the fact that f B u 2 = 1. Thus 



v = f \vu £ \ + [ u 2 <i-c -r 

H^[B^\ Jb 1 JB l 

£' £ £ 



and we have completed the proof. □ 
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We may extend u £ on the entire plane R 2 by setting u £ (r) = for r > K Note that each 
u £ is radially symmetric. Then Lemma |3.8| gives ||w e || H i/ R 2) < K , which implies 



u £ -»■ U weakly in H} (R 2 ) (3.6) 

as e goes to zero (up to a subsequence). Furthermore, by the compact embedding of Hi radial 
functions to Lj functions (cf. Lions paper M), we have 



^UmL 4 r (R 2 ) (3.7) 



as e goes to zero (up to a subsequence). Now, we want to prove that U is nontrivial. By 
Claim Al, we get u 2 £ — In (1 + u 2 ) < \u\. Hence by (|3.6|) , fl3.7|) and Lemma |3l^ (ii), 

— c > liminf /i r e = liminf / |Vu e | + T \u 2 £ — In (l + w 2 )] , (3.8) 

> limm/ / |Vm £ | + \Tu A s > [ \VU\ + \TU\ (3.9) 

e^0+ Jr2 J R 2 

which implies that U is nontrivial. Here we have used again Fatou's Lemma. Otherwise, U = 
and then > — Co > 0, a contradiction. 

To complete the proof, we have to prove that \\U\\ 2 = 1. The idea is to use the concentration 
compactness. We shall prove that the class {u £ } £>0 is neither vanishing nor dichotomy as 
e by contradiction. Suppose {u e } £> o is vanishing i.e. u £ — > in L 2 (R 2 ) as e — > 0. Then 
by Egoroff Theorem, u e — > almost everywhere as e — > (up to a subsequence). On the 
other hand, by ( |3.7| ) and Egoroff Theorem, it E — >■ J7 almost everywhere as e — >• (up to a 
subsequence). Consequently, {/ = but U is nontrivial. This gives a contradiction, so the class 
{u £ } £>0 is not vanishing. 

To show that {u £ } £>0 is not dichotomy as e — > 0, we study the solution profile of u £ . The 
minimizer u £ satisfies the following equation: 

Au £ - r— ^ = A £ n £ in Si (3.10) 

with the zero Dirichlet boundary condition u £ = on 95 1 , where A e is the associated Lagrange 

e 

multiplier. Multiply equation ( |3.10|) by u £ and integrate over B\. Then using integration by 

e 

2 4 

parts and J B u 2 = 1, we get A e = X £ J B ^ u 2 = —J Bl |Vw e | — T J B j^z- Hence by Lemma pT8| , 

£ £ £ £ 

|A £ | < K u (3.11) 

where i^i is a positive constant independent of e. 

Since ti £ = u £ (r) is radically symmetric, equation ( p.lOj ) and the zero Dirichlet boundary 
condition can be reduced to a boundary value problem of an ordinary differential equation as 
follows: 

u" £ + \u' e - = \ £ u £ for < r < i, 

u' £ (0) = 0, u £ (I) = 0. 1 ' J 

From the energy comparison, we may set u £ (r) > for < r < K 
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Lemma 3.9. The minimizer u £ = u £ (r) is positive and monotone decreasing with r. 
Proof. Suppose X £ > 0. Then equation ( |3.10| ) implies 

u 3 

Au £ - \ £ u £ = r— £ — < in Bi 
l + u 2 E ~ £ 

with the zero Dirichlet boundary condition u £ = on dB\. Hence by the strong maximum 
principle, u £ > in Bi and 

e 

u £ (r) = min u £ (x) = min u £ (x) > min u £ (x) = min u £ [x] = u £ (s) 

x£dB r x£B r xGB a x£dB s 

for < r < s < i since u £ = u £ (r) is radically symmetric. This shows that u £ = u £ (r) is 
positive and monotone decreasing with r. Similarly, if A e < 0, then equation ( |3.10| ) implies 

Au £ = (\ £ + r— u £ <0m Bi 

with the zero Dirichlet boundary condition u £ = on dBi. Thus using the strong maximum 

e 

principle again, we may prove that u £ = u £ (r) is positive and monotone decreasing with r, and 
thus complete the proof of Lemma [3.9| . □ 

From Lemma |3l| u £ can not be splitted into two parts as e — > 0+, and hence {u £ } £>0 
can not be dichotomy as e — > 0+- Therefore, by the concentration compactness Theorem, 
u £ — )■ U in L? (M 2 ) as e — » 0+ (up to a subsequence), which implies ||f/|| 2 = 1 since \\u e \\ 2 = 1. 
Now we claim that the limit function U satisfies 

U 3 

Au ~ T T+lP = XoU in r2 ' (3,13) 

U = U (r) > is radially symmetric, and lim U (r) — 0, where Ao is the limit of A £ 's (up to a 
subsequence) since (|3.11|) implies 

A e — > Ao as e — > 0+ (up to a subsequence). (3-14) 

Let (f) e Cq 30 (R 2 ) be any test function. Since u £ satisfies ( |3.1U| ), then 

Vu £ -V<P + T f -^—^=-^1 u £ <\>. (3.15) 

j R 2 i + ui y K 2 



Hence (|3T6|) and ( |3.7| ) give 



Vu £ -V(j)^ / VU-V(f> 

U £ (f) / 



and 



U e i f m f U ^ 
U £ (p 



1 + J R 2 J R 2 1 + U| 
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iR2 1 + Uj Jl2 1+ M| 



Note that 



< 



U + u f 



1 + U 2 

U<f>- 



4 and 

3 



U 



U + u £ 



l + U 2 ) 



Ucf) 



l + U 2 



u 



u 3 



l + U 2 ^' 



U+u e 



< 



(1 + u 2 ) (l + U 2 ) 
Thus and (fTTJ) imply 



(l+tt§)(l+E/ 2 ) 



C/0 



< 



< lic/l 



(i + «2) (i + cr») 



2 since 

1 1 
< - + - 



vt/- v0 + r 



{7 3 



l + U 2 



■Ac 



{7</> 



and then {7 satisfies (|3.13|) and lim {7 (x) = 0. Moreover, U — U (r) > is radially symmetric, 

|a;|— >oo 

since each ii £ is positive and radially symmetric. Therefore, the equation for U can be written 

i u 3 



as 



u" + w - r 



A n {7 for r > 0, 



t/'(0) = 0, C7(oo)=0 . 
By the uniqueness of ordinary differential equations, we may assume U (0) > 0. 

Lemma 3.10. U (r) > /or r > 0. 

Proof. We may prove by contradiction. Suppose there exists r > a minimum point of U such 
that U (ro) = 0. Then U' (ro) = U (ro) = 0. Hence by the uniqueness of ordinary differential 
equations, U = contradicts ||{7|| 2 = 1 and we have completed the proof. □ 



Due to lim U (r) =0, there exists R\ > such that < {7 (r) < 1 for r > R\. By equation 



Aq G 7°°. Hence by the standard 



(I3TT3D , AC/ = (r^ + A ) U G L 4 since 

regularity theorem of Poisson equation, {7 G W 2,4 (B Rl ) and then by the Sobolev embedding 
W 2 ' 4 (B Rl ) C L°° (t3 Ri ), we obtain 



Lemma 3.11. U (r) < K 2 for r > 0, where K 2 is a positive constant. 



Now we prove that A is positive by contradiction. Suppose A < 0. Then equation ( |3.13| ) 
and Lemma |3.10| imply 

U 3 



AU = \ U + T- 



< in 



l + U 2 

Hence by Lemma |3.11| and the Liouville Theorem, U must be a constant function, i.e. {7 = 0, 
which is impossible. Therefore, we conclude that 
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Lemma 3.12. A > 0. 

Since T < 0, the equation ( |3.13| ) becomes 

AU - X n U 



U 3 



1 + U 2 



< in 



Hence, as for the proof of Lemma |3.9| , we may use the strong maximum principle to prove that 
Lemma 3.13. U = U (r) is monotone decreasing with r. 

In this manner, we have completed the proof of Theorem [2.1| (ii), and with this the complete 



proof of Theorem 2.1 



4 Conclusion 

In this paper the existence of ground states in the form of counterpropagating solitons in a 
multidimensional nonlinear Schrodinger model of paraxial beam propagation in media with 
saturable nonlinearity has been proven. From the proof, a threshold condition on the beam 
coupling constant for the existence of such fundamental solitons followed. 
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